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Abstract 

Multifragment disintegrations, measured for central Au + Au collisions at 
E/A = 35MeV, are analyzed with the Statistical Multifragmentation Model. 
Charge distributions, mean fragment energies, and two-fragment correlation 
functions are well reproduced by the statistical breakup of a large, diluted 
and thermalized system slightly above the multifragmentation threshold. 

Hot nuclear systems produced in intermediate energy nucleus-nucleus collisions 
are known to decay by multiple fragment emission [jl], ^, ^. While this decay mode 
may be related to a "liquid-gas" phase transition in finite nuclear systems, both 
statistical and dynamical aspects of multifragmenting finite nuclear systems must 
be understood before inferences about nuclear phase transitions can be made. Since 
the thermodynamic limit should be more readily reached for the heaviest possible 
nuclear systems, studies of central colhsions between heavy nuclei, such as Au + Au, 
are of particular relevance. 

Central collisions between heavy nuclei produce maximum fragment multiplici- 
ties at incident energies of E/A ^ lOOMeV ^, Hot nuclear systems formed at 
this energy have, however, been shown to undergo a rapid collective expansion |^, |^ 
which complicates a statistical interpretation of the decay [H, The dynamics 
of collective expansion should be less important at lower incident energies |p!0| , |11|] 



where comparisons with equihbrium statistical model calculations [0, |^ should be 
more appropriate. In this paper we perform such a comparison for fragments pro- 
duced in central Au + Au collisions at E/A = 35MeV. 

The experiment was performed at the National Superconducting Cyclotron Lab- 
oratory of the Michigan State University. Experimental details have already been 
reported in Refs. [0, |TT[]. Briefly, charged particles of element number Z < 20 were 
detected at 23° < 9iab < 160° by 171 phoswich detector elements of the MSU Mini- 



ball array |]12| ; fragments with charge up to Z = 83 were detected at 3° < Oiat < 23° 
by the Multics array [TBI. The geometric acceptance of the combined array was 
greater than 87% of 47r. The charge identification thresholds in the Miniball were 
Eth/A ^ 2,3,4Me\/ for Z = 3,10,18, respectively, and Eth/A ^ l.bMeV in the 
Multics array, independent of the fragment charge. For the present analysis, central 
collisions are selected by requiring observed charged particle multiplicities Nc > 24 
(representing about 10% of the total reaction cross section |1TT|). 

For central events, the fragment emission was found ^ compatible with 
a near-isotropic decay of a source consisting of more than 300 nucleons and with 
negligible contributions from the decay of projectile and target-like residues. Here 
we address the question to which degree fragment emission can be described by a 
statistical equilibrium calculation performed with the Statistical Multifragmentation 
Model (SMM) of Refs. |P, This model was successfully applied at higher energies 



to the interpretation of Aw-induced projectile fragmentation reactions for which 
collective expansion velocity components have been shown to be negligible and also 
for central collisions for which the model has to be modified to take into account 
the large observed expansion velocities. 

The model is based upon the assumption of statistical equilibrium at a low- 
density freeze-out stage of the reaction at which the primary fragments are formed 
according to their equilibrium partitions. The equilibrium partitions are calculated 
according to the micro-canonical ensemble of all break-up channels composed of 
nucleons and excited fragments of different masses [|l6l. The model conserves total 
excitation energy, momentum, mass and charge number. The statistical weight 
of decay channel j is given by Wj oc exp Sj{E*, A^, Zg), where Sj is the entropy 
of the system in channel j and E*,As, and Zg are the excitation energy, mass 
and charge number of the source. Different breakup configurations are initialized 
according to their statistical weights. The fragments are then propagated in their 
mutual Coulomb field and allowed to undergo secondary decays. Light fragments 
with mass number Af < 4 are considered as stable particles ("nuclear gas") with 
only translational degrees of freedom; fragments with Af > 4 are treated as heated 
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nuclear liquid drops. The secondary decay of large fragments {Af > 16) is calculated 
from an evaporation-fission model, and that of smaller fragments from a Fermi break- 
up model The simulated events were then filtered with the acceptance of the 
experimental apparatus |[TT| , p!7| . The same normalization to the total number of 
events was applied to the experimental and calculated distributions, which thus 
may be compared on an absolute scale. 



In its original version [Q, the SMM only incorporates thermal degrees of free- 
dom, i.e. the fragment energy distribution was determined from the source temper- 
ature and then modified by final-state Coulomb interactions and secondary decays. 
In the present work, we also allow for a collective radial expansion of the system 
which could arise from a rapid thermal expansion, possibly aided by an initial com- 
pression. Specifically, we assume that modest collective velocity components do not 
influence the fragment formation probabilities for a given thermal energy; this as- 
sumption is reasonable for E flow I A < 3MeV A self-similar collective expansion 
was assumed, vjiow oc r, where r is the distance from the source's center of mass. 
This collective velocity was added to the thermal fragment velocity. The energy 
balance was taken into account. 

In order to obtain an equilibrium freeze-out condition, we need to estimate mass, 
charge, and energy carried away by particles of early emission. For this purpose, we 
assume that the preequilibrium emission consists primarily of n,p,d,t,^ He and a- 
particles distributed uniformly in the available phase space in centre-of-mass system. 
Following Ref. 0, we search for parameters (mass Ag, charge Zg, excitation energy 
E*, and freeze-out density ps, assuming Zg/Ag = 79/197) of thermalized sources 
which can fit the experimental data. A rough estimate of Ag, Zg, and E* was 
obtained by comparing measured and predicted multiplicities of charged particles, 
Nc, and intermediate mass fragments Njmf ^ Z < 30). To explore sensitivities 
to the assumed freeze-out density, we performed calculations for pg = po/S and 
Po/6 (where po ~ 0.15/m~^ is the density of normal nuclear matter). Assuming 
Ps = Po/3, the model can reproduce central collision data for Zg ^ (0.7 — 0.9)Ztot 
and E* ^ (0.6 — 0.7)Etot where Ztot and Efot denote the total available charge 
and center-of-mass energy. For pg = po/6, the source parameters are 10% smaller. 
The approximate temperature and entropy per nucleon of the extracted source are 
T ^ 6MeV and S/A ^ 1.5 - 1.6. 

Improved source parameters can be obtained by analyzing the inclusive charge 
distributions N{Z) = Yield{Z)/Nevents (Fig- 1 a,b) and the exclusive charge distri- 
butions for the "first" six fragments ordered according to Zi > Zk if i < k (Fig. 2). 
We show SMM calculations for two sets of source parameters: 
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As = 343, Zs = 138, E*/A = G.OMeV, = po/3, 

A, = 315, Zs = 126, E*s/A = A.SMeV, = po/6. 

In the second case an additional radial flow energy Ejiow/A = O.SMeV was intro- 
duced (the total energy of the source was thus Es/A = 5.6MeV). This modest 
amount of radial flow is in agreement with the extrapolation of data at higher ener- 



gies to E/A = 35MeV [|1^]. As shown below, both sets of calculations well reproduce 
the measured observables. We have also considered a collective rotational motion 
with the same energy per nucleon as the flow energy and we found our main conclu- 
sions to be essentially unchanged as compared to the case of collective expansion. 

In Fig. 1, the sensitivity to the excitation energy is illustrated by additional cal- 
culations performed for E*/A ± IMeV per nucleon (dot-dashed and dotted curves). 
A comparison of flltered and unflltered calculations (solid and dashed curves, re- 
spectively) shows that distortions of the Z-distribution from to the experimental 
apparatus are relatively small, except in the region Z > 20 where the Miniball de- 
tectors lose charge resolution. In the region of large Z, the charge distribution falls off 
more steeply than expected for exponential or power-law distributions. This steep 
fall-off, reproduced in the calculations, is an effect of charge (mass) conservation for 
a flnite system. 

One can readily understand why sources of different excitation energy and density 
can produce similar fragment distributions. A system of lower density has smaller 
Coulomb barriers for fragment formation and thus requires a smaller temperature 
to break up into fragments. During the expansion to a lower density additional 
particles will be lost and some of the internal energy may be converted into radial 
flow. Therefore, the low-density source should have less mass and excitation energy. 

Mean values, {E/A) (solid points), and standard deviations ^E/A (vertical bars) 
of the kinetic energies per nucleon of fragments emitted at 6cm = 90"± 10° are shown 
in Fig. 3 (left panel) as a function of Z. The solid and dashed curves show the results 
of SMM calculations at {E/A) ±AE/A for ps = po/3 and po/6, respectively, flltered 
by the acceptance of the experimental apparatus. Right panels of the flgure show 
the dependence of {E/A) on O^m for Z = 6, 10, and 14. (The small rise in {E/A) at 
forward angles and the small dip at backward angles are caused by the acceptance 
of the experimental apparatus [|ll|]-) Overall, both SMM calculations reproduce the 
data rather well - the smaller Coulomb repulsion from the lower-density source is 
compensated by the added collective expansion energy. Both calculations underpre- 
dict the kinetic energies of lighter fragments {Z < 5) and of fragments emitted at 
forward angles, possibly due to the presence of some nonequilibrium fragment emis- 
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sion in the data. The bulk of the data is, however, consistent with near-equihbrium 
emission from a single source |Tl|| . 

In order to test the predicted spatial separation of emitted fragments, we have 
constructed two- fragment correlation functions ||T9|, 



^back (^red) 

where 



^3 



is the "reduced" relative velocity of fragments i and j (i 7^ j) with charges Zj and Zj] 
Yivred) and Yi,ack{vred) are the coincidence and background yields for fragment pairs 
of reduced velocity Vred] and C = Nback/Ncoinc where Ncoinc and Nback are the total 
number of coincidence and background pairs. The background yield was constructed 
by means of the mixed event techniques pO|, |2l| . 

Because of the large differences in dynamic range and resolution between the 
Multics and Miniball arrays, we separately evaluated the two-fragment correlation 
functions measured with the two devices. Figures 4 a) and b) show, respectively, 
two-fragment correlations constructed solely from fragments detected in the Mul- 
tics array {3 < Z < 30, 8° < Oiab < 23°) and in the Miniball {3 < Z < 10, 
23° < Blab < 40°). The solid and dashed curves show the results of SMM calcu- 
lations performed for ps = po/S and po/6, respectively, filtered for the acceptance 
of the experimental apparatus. Overall, the measured correlation functions are ap- 
proximately reproduced by the calculations. Some discrepancies may be caused 
by the small differences found in the event charge partitions shown in Fig. 2 (see 
analysis on Ref. ||22[). The increased overshoot just after the Coulomb hole for the 
case Ps = Po/6 is caused by the flow: more "organized" motion produces additional 
correlations in single events. 

Two-fragment correlation functions measured in this experiment are also re- 
produced fill by assuming sequential emission from the surface of a spherical 
source of charge Zg = 138, breakup density ps = Po/4, radial expansion velocity 
Vfiow = 1-4 cm/ns and average interfragment emission time r ^ 85/m/c. This 
ambiguity in interpretation is due to the fact that the angle integrated correlation 
functions are affected by both the source's size and its lifetime, resulting in the well 
known space-time ambiguity [0, |2^. However the energy and angle dependent anal- 
ysis ||2^, |2^, ^ indicates that fragment emission may be sufficiently fast to justify 



the approximation of instantaneous emission. In the SMM we increase effectively 
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the time of fragment production (and decrease the initial correlation) by including 
the secondary de-excitation. 

The source parameters extracted in our analysis are consistent with expectation 
from dynamical simulations. We have performed such simulations with the Boltz- 
mann Nordheim Vlasov (BNV) model |^ for central collisions {b < 1 fm) using 
the mean-field approximation with a soft equation of state. Approximately 100/m/c 
after the initial contact between projectile and target, the calculations predict the 
formation of a single source of mass A ^ 324, charge Z ^ 136, density p = po/2, ex- 
citation energy per nucleon E* /A ^ 6MeV, and rotational energy Ej-ot ~ lOMeV. 
Unfortunately, BNV-type models do not treat emission of single fragments. 

In conclusion, fragment emission observed in central Au + Au collisions at 
E/A = 35MeV is largely consistent with the statistical break-up of a single source 
of excitation energy per nucleon E*/A ^ 5 -r- 6MeV and density ps = po/3 -ir po/6. 
Collective expansion energy is not excluded but it is smaller than ~ IMeV/u. 
These conditions are comparable to those observed for peripheral projectile frag- 
mentation reactions at higher energies {E/A > AOOMeV). However, these reactions 
can only produce smaller sources of mass and charge numbers Ag < Aprojecuie and 

Zg ^ Z projectile- 
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Figure captions 

Fig.l: Charge distribution N{Z). Points show experimental data and Unes show results 
of SMM predictions for sources with parameters Ag = 343, Zg = 138, /A = 6.0MeF, 
= pq/3 (part a)) and A^ = 315, Zg = 126, E*/A = A.SMeV, EjioJ/A = O.SMeF, 
Ps = Po/6 (part b)). Dashed curves are the unfiltered calculations and solid curves are the 
filtered ones. The dot-dashed and dotted curves represent filtered calculations for thermal 
excitations E*/A + IMeV/u and E*/A — IMeV/u, respectively. 

Fig. 2: Charge distributions of the six heaviest fragments in each event, ordered such as 
> Zy; \i i < k. Experimental data are shown by points, the solid and dashed curves 
show the results of SMM calculations for ps = po/S, and ps = po/Q, respectively (other 
source parameters as in Fig. 1). 

Fig. 3: Mean centre-of-mass kinetic energy per nucleon, {E/A), as a function of the 

charge Z, for fragments emitted at 9cm = 90° ± 10° (left panel) and (for Z = 6, 10, 14) as a 
function of 9cm (right panels). Points give the experimental values of {E/A) and vertical 
bars give the standard deviations ^E /A of the distributions. The solid and dashed lines 
are SMM predictions of {E/A) (in the left panel ^IS.E/A) for ps = po/S, and ps = po/Q, 
respectively (other source parameters as in Fig. 1). The energy range is the same in the 
left and in each right panel. 

Fig. 4: Two-fragment correlation functions 1 + R{vj.ed) for 3 < Z < 30 and 8" < Oiab < 23° 
(part a)) and for 3 < .Z^ < 10 and 23° < 9iab < 40° (part b)). Full points show experimental 
data. The solid and dashed lines are SMM predictions for ps = po/3, and ps = po/6, (other 
source parameters as in Fig. 1). 
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